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Abstract 



An active elastodynamic cloak destructively interferes with an incident time harmonic 
in-plane (coupled compressional/shear) elastic wave to produce zero total elastic field 
over a finite spatial region. A method is described which explicitly predicts the source 
amplitudes of the active field. For a given number of sources and their positions in 
two dimensions it is shown that the multipolc amplitudes can be expressed as infinite 
sums of the coefficients of the incident wave decomposed into regular Bessel functions. 
Importantly, the active field generated by the sources vanishes in the far-field. In practice 
the infinite summations are clearly required to be truncated and the accuracy of cloaking 
is studied when the truncation parameter is modified. 
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^ ■ 1 Introduction 

m 

If} ■ The main function of a cloaking device is to render an object invisible to some incident wave 

as seen by some external observer. Over the past decade, a great deal of effort has been 
focused on passive cloaking, using metamaterials to guide waves around specific regions of 
space, see e.g. the highly cited works [H [21 [3]. In recent times a rather different approach 
to cloaking has been noted as an alternative. It has been named active exterior cloaking 
and it relies on a set of discrete active sources, lying outside the cloaking region, to nullify 
the incident wave whilst their own radiated field must be negligible in the far-field. Interest 
has focused on the Helmholtz equation in two dimensions [H [5j El [HJ [9] . In the work of 
Vasquez et al. [6] Green's formula and addition theorems for Bessel functions were used to 
formulate an integral equation, which was then converted to a linear system of equations for 
the unknown amplitudes. Crucially, the integral equation provides the source amplitudes as 
linear functions of the incident wave field. It was shown that active cloaking can be realized 
using as few as three active sources in 2D. Further work to render the linear relation for 
the source amplitudes in more explicit form was developed in [8] and extended to the three 
dimensional Helmholtz case in [7]. In jS], the integral representations of Vasquez et al. [8] for 
the source amplitudes were reduced to closed- form explicit formulas. This obviated the need 
to reduce the integral equation of Vasquez et al. [6] to a system of linear equations which 
are then required to be solved numerically or to evaluate line integrals, as proposed in [8]. 

There is, of course, a strong link between active exterior cloaking and the notion of anti- 
sound or in the context of elastic media anti-vibration. Interestingly the notion of anti-sound 
appears to have been considered first in a patent published in 1936 by Paul Lueg, see e.g. [10] ■ 



* Corresponding author: William.Parnell@manchester.ac.uk 



1 



Active elastodynamic cloaking 



November 26, 2012 



2 



The subject has focused greatly on the desire to reduce the magnitude of a radiating field or to 
create so-called quiet zones in enclosed domains such as aircraft cabins using simple sources. 
The idea to suppress completely the sound field in a finite volume inside an unbounded 
domain using the Kirchhoff-Helmholtz integral formula and thus employing a continuous 
distribution of monopoles and dipoles is described in [11]. Anti- vibration techniques have 
also been developed [121 03] • I n general the focus of anti-sound is to reduce the sound 
radiated from a sound source or to create a zone of silence by employing a finite number of 
radiating sources. The active field is not required to be non-radiating however. Furthermore 
very little work in the anti-sound community has focused on the exact shape of the quiet zone 
with the exception of |14| who calculated, numerically the zone of silence (< lOdB) region 
created when the amplitude of a single secondary source was chosen to reduce noise due to a 
single primary source. 

The aim of active exterior cloaking is to render the total field zero inside some prescribed 
domain (the cloak or zone of silence), whilst ensuring that the active field itself is non- 
radiating. The technique introduced in the early active exterior cloaking work enables a 
cloaked region to be identified clearly by the use of Graf's addition theorem. This approach 
allows precise determination of the necessary source amplitudes. 

The infinite series associated with the multipole expansion of the rath active source is 
formally divergent inside the circle that is centered on source itself, i.e. for |x — x m | < a m 
in the notation used later on. Therefore the representation for the source is not valid in the 
domain in which it resides! This point has not been stressed in the active cloaking community, 
although a related point was noted in the anti-sound community in [15]. Clearly this point 
motivates the truncation of the series which has to be done practically in any case. This 
limits accuracy but as we shall see later in many cases, only a small number of multipoles 
are required. 

As yet it does not appear that active exterior cloaking has been applied to the elastody- 
namic context. This paper will focus on the relevant two dimensional active elastodynamic 
cloaking problem. In general, elastodynamic cloaking problems are more difficult to study 
than their acoustic or electromagnetic counterparts. Indeed in the case of passive elastody- 
namic cloaking, this is due to the lack of invariance of Navier's equations under coordinate 
transformations [16] unless we relax the minor symmetry property of the required elastic 
modulus tensor. The latter can be achieved by using Cosserat materials |17[ [18] or by em- 
ploying nonlinear pre-stress of hyperelastic materials |19[ [20l 121 j . Here we show how the 
active approach to cloaking can be employed in the elastodynamic case for the fully coupled 
two-dimensional (in-plane) compressional/shear wave problem. As in the approach of [9] 
we write down the relevant integral equation by employing the isotropic Green's tensor in 
this case. The required source amplitudes for arbitrary wave incidence can be determined 
explicitly by using Graf's addition theorem. 

We shall begin in Sj2] with a statement of the problem, a review of the governing equations, 
and a summary of the main results. The relevant integral relation is derived in from which 
the main results regarding the explicit form of the source amplitudes are shown to follow. 
We consider both compressional and transverse (shear) wave incidence. We also describe the 
form of the active source field and the issues associated with divergence described above. 
Numerical results follow in SJ5J 

2 Problem formulation and main results 
2.1 Problem overview 

Let us consider the two-dimensional configuration where the active cloaking devices consist 
of arrays of point multipole sources located at positions x m 6 M 2 , m = 1,M as depicted in 
Figure [TJ These sources can give rise to both shear and compressional elastic waves. The 
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active sources lie in the exterior region with respect to the cloaked region C and for this 
reason, this type of cloaking is called active exterior cloaking [5]. Objects are undetectable 
in the cloaked region by virtue of the destructive interference of the sources and the incident 
field with the result that the total wave amplitude vanishes in the cloaked region C. As 
described in [9] this gives rise to three significant advantages over passive cloaking: (i) the 
cloaked region is not completely surrounded by a single cloaking device; (ii) only a small 
number of active sources are needed; (hi) the procedure works for broadband input sources. 
The principal disadvantage of the method is of course that the incident field must be known. 

The M active sources give rise to a cloaked zone C is indicated in Figure Q] by the shaded 
region whose boundary C is the closed concave union of the circular arcs dC m m = 1,M, 
{a m , #[ m \ flij } associated with the source at x m . In the general case {a m , , 0^} are 
distinct for different values of m. Note that the wave incidence shown in Figure [T] is a plane 
wave although the solution derived below is for arbitrary incidence. We therefore have to 
determine the amplitudes of the active sources as a function of the incident wave, and then 
prove that the cloaked region is indeed the closed region C as indicated in Figure CD Let us 
also define the notation A m as the circular domain of radius a m that contains the mth active 
source at its centre. We also define the union of these domains A = L)5f =1 ./4 m . 



Figure 1: Insonification of the actively cloaked region C generated by M active point multipole 
sources at x m , and active sources regions A m , m = 1,M. The incident field in this case is a 
plane wave with wave vector k in the direction ip. 

2.2 P/SV in-plane wave propagation 

We consider time harmonic solutions with the factor e~ lLut understood but omitted. Navier's 
equations in two dimensions for the displacement u = (ui,U2), Uj = Uj (xi,X2), are, 



where f is the forcing, Cijki = \5ij5ki + fJ-(o~ikO~ji + b~u5jk) and the summation convention on 
repeated indices is understood. Hence, in the absence of forcing (f = 0), 




djO-ij + puj Ui = fi, aij = C^ki diu k , 



(1) 



(A + 2jj) VV • u + fiV 2 u + pu 2 u = 0. 



(2) 
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The Helmholtz decomposition for the displacement, 

u = V$ + Vx(^k) (3) 

leads to separate Helmholtz equations for the scalar potentials 

V 2 $ + kl§ = 0, V 2 ^ + k 2 s ^ = (4) 



where k p , k s are the longitudinal and shear wavenumbers, respectively: k 2 = uj 2 p/{\ + 2p), 
k 2 = uj 2 p/fi. We also define, for later use, k = k s /k p , or equivalently k 2 = 2(1 — u)/(l — 2v) 
where v is Poisson's ratio. We seek the total wave field in the form of an incident wave, Uj, 
plus the active source field, u^, such that 

u(x) = Ui + u d $(x) = $i + $ d , ¥(x) = + (5) 

We assume the general form of an incident field in the regular basis, and hence 




where the functions U^{z) and V^{z) are defined by 

U±(z) = J n (|z|)e ±i '"*«», V±{z) = H^(\z\)e ±m ^. (7) 

Here argz € [0, 2ir) and arg(— z) = argz ± tt € [0, 2ir). Define the derivative functions 
(z) and V^ 1 (z) as 

I7±'(z) = J' n (\z\)e ±in ^, V±'{z) = H<p\\*\)e*™*\ (8) 
Note that the functions U^{z) and V r ^ 1 (z) possess the properties 

U±(-z) = (-l) n U±(z), V±(-z) = (-l) n V±(z). (9) 
In the following we write Uq and Vo, with obvious meaning. 

2.3 Summary of the main results 

Here we shall state the main results and the required source amplitudes to enable perfect 
active cloaking together with necessary and sufficient conditions on these amplitudes. The 
latter ensures we can compare accuracy of the cloaking technique. We shall prove these 
results in $3j Let {o m , #J m \ 6^} define the circular arc dC m of the closed boundary of the 
cloaked region associated with the source at x m . The active source amplitude coefficients for 
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Derivation of equation (|10p is given in Section 13.41 Alternatively defining a vector a^ = 

a m e(#j m ^)(i = 1, 2), and incorporating equations (|7|) and flS}, equation (jlObp reduces to the 
form 



' ' B m,ln\ i 00 f 

„ =^2 E (-1)«|C^(A^ 

\ m, In/ 



V(a,/3) 



Fi(A; p a, k s a) 
V 2 (k p a, k s a) 



-V 2 (k s a, k p a) 
Vi(fc s a, k p a) 



A™) 

Amy 



, where 
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(m) (m) TT -t, 



U~(a) 



(11a) 
(lib) 



Next, we note that the active source coefficients B^ J , and B^ t must satisfy the necessary 
and sufficient conditions to ensure active cloaking inside the domain C. 

'^X-z(VO =0, 



M oo 



\/n € Z : ^ ^ x< 



m=l i=— oo 



B m,l V n-l(kpXr, 



0. 



(12) 



The first pair of conditions are required to ensure zero radiated field outside the union of 
the active regions, A, and the second pair ensures that the total field is zero inside C. These 

(p) (s) 

constraints on Bz/, and Bi/, will be used to estimate the error in the active cloaking region 
in the following sections by truncating the infinite sums in (|12p . 
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3 Derivation of the source amplitude expressions and con- 
straints 

Let us first formulate the problem in terms of an integral equation. 
3.1 Green's tensor and integral equation formulation 

Consider the particular solution of Navier's equations ([I]) in the presence of a point force, 

f = F5(x-x ) =>- Ui = G %k F k (13) 
where Gi k is the two dimensional (in-plane) Greens tensor. Specifically, Gj&(x) satisfies 

^ijkj + P^Gik = <5jfc5(x), Sjjfc = CijpqGpk } q, (14) 



with solution 



where 



G ik = ipu 2 ) 1 



5 ik k 2 s G s + did k (G s -G p ) (15) 



G s = ^Vb(£ s x), G p = ^V (k p x). (16) 

The solution (|15j) can be checked by substitution into the governing equation (|14p and using 
the identities (V 2 + k^)G a = <5(x), a = p,s. 

It is convenient to work without subscripts, writing (|15p as 



-pw 2 G(x) = VVG P + (IV 2 - VV)G S = VVG p + (Vxk)(Vxk)G s for x / 0, (17) 

where (Vxk)j = e^dj. Using (|13p in the form u = G • F, combined with the Helmholtz 
decomposition for u gives 

V$ + (Vxk)^ = [pu?)~ l [VVGp • F + (Vxk)(Vxk)G s • F] , x + 0, (18) 

implying 



-pu 1 ^ = F • VG P = F&Gp + F 2 d 2 G p , 
-puH = F • (Vxk)G s = F^Gs - F 2 d 1 G 



x + 0. (19) 



This makes it clear that for a standard point source, regardless of the choice of F, both 
compressional and shear waves propagate away from the point source. 

With knowledge of the Green's tensor we can now develop an integral equation for the 
displacement. Indeed, if u is a solution of the homogeneous equations in an infinite domain 
containing a finite region D and a is the associated stress, then by definition of the Green's 
tensor, 

I d5 , n i [uj(y)S ijjfc (y - x) - a i:j (y)G jk (j - x)] = \ ™ k *^ * £ ^' (20) 

J 3D 10, KfD. 



Equation (I20p holds for both Uj and separately inside the cloaked region, since both are 
assumed to be regular there (this is a definition of exterior cloaking). Also, by its definition 
the total field is zero inside the cloaked region with boundary dC, and therefore 

u d (x) = -/ dSn- [u i (y).E(y-x)-o- i (y)-G(y-x)] > x e C. (21) 

JdC 

This is the fundamental relation used to find the source amplitudes. 
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3.2 General expressions for the source amplitudes 

Following the procedure for the Helmholtz problem [9], we first substitute the assumed form 
of Urf into the left member of (|2ip. Then we partition the integral in the right member into 
M segments over {dC m ,m = 1,M} and identify each line integral with the m th component 
of Ud, i.e. the part of the source field from the multipoles at x m . Thus, 

AI f 

= E / dSn ■ My) • s (y - x ) - ^fr) ■ G ^ - x )] 

m=l ^ J dCm 

+ f; te) n VC(fc p (x-x m )) + BW„Vxky/(fc s (x-x m ))H ) xeC. (22) 
We now use the generalized Graf addition theorem [22j eq. (9.1.79)], 



.+ (y _ x)= f (K + (y)^(x), |y|>|x|, 



(23) 



The idea is to write S(y — x) and G(y — x) in (f22j) in terms of sources at x m . This suggests 
using (j23|) for y — x = (y — x m ) — (x — x m ) subject to |y — x m | < |x — x m |. Hence, using 

OH), 



G(y-x)=^ J] {vVf/-(fc p (y-x m ))K + (A ;p (x-x m )) 

n=— oo 

+ (Vx k) (Vx k)U~ (k s (y - x m )) V;+ (fc,(x - x m )) J . (24) 

By virtue of the dependence of the Green's function on y — x the derivatives VV can be 
understood as VyV^ or V X V X or — V^V^., with the same equivalence for (Vxk)(Vxk). 
Inspection of ([22]) suggests that the forms —V y V x and — (V y x k)(Va;X k) are appropriate. 
Taking into account the negative sign in VV — > — VyV x , the Green's function can be written 
in the form 

OO f 

G(y " x) = 4^7 E \ v v u n{k P (y - x m ))v :c y n + (/c p (x - X m )) 

" n=— oo 

+ (V y x k)U~ (k s (y - x m )) (V x x k)V n + (k a (x - x m )) J . (25) 

Substituting from ()25[) into (|22p . and identifying the coefficients of VV^ + (A; p (x — x m )) and 
Vxk^ n + (fe(x - x m )), yields 

B $n = jte\ d Sn • My) • w n (y - x -)) - u *(y) ■ a(p) (My - x -))] . ( 26a ) 

B$n = ^ I d5n- [<7,(y) • (Vxk)[/;(fc s (y - x m )) - u 8 (y) • ff W(fc s (y - x m ))] , 
Apu 2 J aCm 

(26b) 

where 

a ij (^p(y ~~ x m)) = Cij pq U npq (k p {y — x m )j , 

(fc 8 (y - x m )) = C ijpq e pr3 U~ rq (k s (y - x m )). 

Therefore, given the incident field, we are now able to evaluate the required source am- 
plitudes that guarantee zero total field inside the domain C. We can however, make further 
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progress on the integrals in (|26|) in order to render them in simpler form, by using the fact 
that dC m is the arc of the circle of radius a m centered at x m , which is the origin of the shifted 
coordinates y — x m . The integration is therefore simplified using polar coordinates centered 
at x m , combined with the expressions for the displacements and traction components in polar 
coordinates given in terms of the potentials, 



a rr = -\k 2 p <& + 2//($ rr + rfl - -^,b) > 



(27) 



The four distinct terms in the integrals of (j26H . such as d Sri • <Xi(y) • VC/ n (k p (y — x m )) , then 
follow by identifying $ — > U~ (k p a), ^ — > U~ (k s a), where a(9) = y — x m is the radial vector 
of constant magnitude a m . Thus, 



d Sn ■ <r,j • VU n = d 5 



d 1 d 

a ^rjrU~ (fcpa) + <Tir6--^U n (k p &) 



dr n V P ; ' uo rd9 
d9 (oir r kpa m U~' \k p a) - ino ir eU~ (k p a 
(p)i 

r6 J 



(28a) 



d Sn • <j( p ) • u.j = d S \ui r o)P) + UiQcr, 



dSn-ov (Vxk)£/~ 
d 5n • er^ • u, 



d 9-^— | 



dS 



+ ti^ 2in[U~(k p a) - k p a m U n ~' (k p a)]j , 

1 9 9 
a irr-ggU n (k s a) - o ir6 —U n (k s a) 

-d9^ina irr U~(k s a) + /c s a m cr.j r 6iC/ r 7'(A; s a)^ , 

d5[ui r cr^ + n ie c7^] 

dO— (u ir 2in\U~(k s a) - k s a m U~' (k s a)] 
a m V 

+ ^[((fcsam.) 2 - 2n 2 )U-(k s a) + 2/c s a fn [/- / (fc s a)) 



(28b) 



(28c) 



(28d) 



Noting the reversal of the sense of the integral in equation (|26p and incorporating equation 
5aJ) leads to 
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+ 2/ 



j ; (4 m) )-4 m) j/(4 m )) 
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,(m) 



d0e 



iai m )j/(a( m ))^ + ZJK«i m) ) 



(X; 
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2/ 2 )J;( 



a 



(m) N 



/' 

+ 2ai m) J/(ai m ); 



u-ie 



+ 21 



(29a) 



u 

a. 

(29b) 



where a 
„( m ) 



(m) 



k p a rn , = k s a m , 9\" 1 ' and O^""' are the angular positions of the vectors 

i = 1,2, which describe the initial and final positions of segment dC m . 



,(m) 



,(m) 



a m e(9, 
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Equations (|29p provide expressions for the source amplitudes for any time harmonic incident 
field. 

Let us now specialize the result to the specific case of plane wave incidence. This is 
important in its own right but also allows us to derive the general incident wave case by 
integration as we shall show. 

3.3 Plane wave incidence 

Let us define 

u^(x) = e ifc ^)- x , a = p,s (30) 

where e(ip a ) = (cosiJj a ,smiJj a ) so that u^ a correspond to compressional (p) and shear (s) 
plane waves of unit amplitude. 



3.3.1 Longitudinal incident plane wave 

Consider now longitudinal plane wave incidence 

$i(x) = A p u^ p (x) 



(31) 



where A p = const is a known wave amplitude. Then using the relation ^(y) = <&i(x m )it,/, (a) 
with a = a m e(9), and eq. (|27j) with $ = <E»j, VI/ = 0, reduces equation (I29p to the form: 



B 



( p ) _ ^(X,; 



m, I 



AkoIs^ Jo™ 



,(m) 
'2 



(32a) 



m, £ 



+ /4 m ) J;(4 m) ) sm2(9 - if; p ) - i2lsm(9 - i/j p ) \ji(a^) - J{ '(a p m) ) 
- cos(9 - W) [(«i m)2 - 2/ 2 ) J/(4 m )) + 2«M J/(4 m) )] |, 

! dfle"*' V( a ) • | t 2sin2 ( " ^) " Z<4 m) JK«i m) ) 
zsin2(0 - ^ p )a(, m )a( m ) Ji'(ai m) ) + i2 Z cos(0 - Ji{a^) - a^J/^) 
sin(0 - Vv) \{a ( r )2 ~ 2l 2 )j l {a { [ n) ) + 2a^ J/ (a^ )]} ■ (32b) 



Then noting that (a) = e *fcpOm «>s(0-v<p) _ e i<4 m) cos(0-^ p ) ; e q Ua tion ([32]) can be written 



(p) _ z$i(x m ) 



m, I 



Ana 



(m) 



• I af) 2 J/(a p m )) [2L' '(4 m )) - {k 2 - 2)L (a^) 



-2la^J l (a^)L' 1 (a^) - 2lL 1 {c$ £ >) [ji(a^) - a^J^a^) 
+ L '( a ^) \{at )2 ~ 2l 2 )j l {a^) + 2a^ Ji'(a^) 



(33a) 



B 



( s ) _ $j(x r 



m, i 



4Ka 



(m) 



^ - /af) J ; («i m) ) [2L "(4 m )) - {k 2 - 2)L (aH) 



2 «M J/ (ai m ) )L; (qM ) + 2 i L{, (a p m ) ) [ J, (a( m ) ) - qH j{ ) 
L!(a p m )) [(a^ 2 - 2/ 2 ) J z (ai m )) + 2aHji'(aH) 



(33b) 
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Lj(a) 



l Vv 



dfl(sin#) J e 



j p i(a cos 0—10) 



, J = 0,1. 



(34) 



L (a) can be evaluated by using the Jacobi- Anger identity e txsm& = Yl^=-oo Jn{x)e m9 , 



Lo(a)= V J„(a)i" / ^ dfle"*^ 



n=— oo 
oo 



„i(n+l)ib v /•„% 



(m). 



(35) 



Integration by parts yields -Li(a) in the form 



T /s t r„\ ^ i(acos0-W) 

a ia 



(36) 



Taking into account the Jacobi- Anger identity and equation (|35p . the function L\{a) and its 
derivative L[(a) can be expressed 



Li(a) = - y~] J n ( 



aim 



n=— oo 
oo 



n + l 



(37a) 



L i(*) = \ E -* n+ >^(«) " '»(«)] ^— 7- • [e-^+<)4 m) _ ^("+0^]. (37b) 
or ' L J n + 1 1 J 



Introducing the explicit results for the functions Lq{q) and L\{a) into (|33|) yields expressions 
for the amplitude coefficients in the form: 
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(m) " l ^P 
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2J/'(4 m) ) 



+ 



2/^ 



a 



(m) 



Ji(4 m) ) 



: V«(^ m) )[4ffJi(4 TO) )- J i / ( t 4 ro) ) 



n 



(m) ' 



a 



(38a) 



(s) _ $ t (x m ) ^ jg+i e i^ P J / 



E 



4k ^— ' o + Z k 

g=— oo v 



JK«i m) ) 



-2gJ/(a( m )) 
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2j/(aM)-(K 2 -2)J 9 («f)) 
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Jq(al )) 
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al m) -2/ 2 



JK4 m) ) + 2J/(ar ; ) 



a: 



(m) 



. [ e -i(q+l)e^ _ e -i(q+l)0[^y 



(38b) 
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After some simplification eq. (|38p can be written as 




E 



9»<#p . \ P -i(g+i)e ( 2 m) _ p -i(q+i)e[ m} - 



4a s q=-oo 



v -;[ a M 2 -2Z?]J,(ai m) ) -24 m) ai m) J/(ai m) ) y 



/^(4 m) A 



(39) 



3.3.2 Transverse plane wave incidence 

Consider now an incident transverse plane wave 

9i = A, e tt '*(iM-* 



(40) 



where A s = const is a known transverse wave amplitude. Entirely analogous calculations to 
the compressional wave case yield the source amplitudes in the form 



TO, I 

V -7? (p) 



(-"-mj 



4a 



(m)' 



E 

c/= — oc 



(m). 



\ —i[a 



(m.) (m.) ,/, (m)^ 



(41) 



3.3.3 Plane wave incidence summarized 

Adding the separate results of eqs. (|39p and (|4ip gives for combined incidence 
the source amplitudes 



(42) 




'ml 2 Z-/ L 

(m) (m) 



(m). 



/I ( m ) 

4«s <?=-oo 



(/ I'm i 
«i(a P ,a£ ) 
/ (m) (m)x 
■u 2 (ap ') 



+ * i (x m )e*^« 

where the vector v(a,/3) = (ui, i^) 7 " is defined in (|10cj) . 



—V2(a s ,a p ) 

/ (m) (m)\ 

vi(al ,a p ') 



(43) 



3.4 Arbitrary incident field as superposition of plane incident waves 

The general form of incident field given by equation (|6ap can be constructed as superposition 
of plane incident waves of the form (jUJ). This will enable us to find the general form of 
the amplitude coefficients for incident waves of general form as a superposition of solutions 
for plane waves given by (|43p . Recall the incident field for a combined incident plane wave 
having the form: 



^(x) 



p ik s e(4>s)-> 



E 

q= — oo 



iie-^U + ikpX) 
iie~ iq ^U+(k s x) 



(44) 
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Multiplying the first row of equation (jUJ) by (z~( n+ '?) /2ir)e l ( n+q ^ p and the second row by 
(i~( n+ i) j2'K)e L ( n+q ^ s , integrating with respect to ip p and ip s respectively between and 2tt 
and then evaluating at x = x m we find 



l~{n+q) 1-2-K 







2tt 



2tt 







d^^x™)^™^ = U+ +q (k p x m ), 



(45a) 
(45b) 



To obtain the form of the amplitude coefficients given by eq. (|10p for the general incidence 
(f6a|) we multiply the first and second of equations (03]) by An^ n (i/j p , ip a ) and r y n (ippj ips) 
respectively, where jniipp^s) = i _2n /(27r) 2 e m ^ p e m ^ s , carry out the double integration with 
respect to ip p and ip s between and 2tt, incorporate (j45|) and sum over all n £ Z. 

3.5 Necessary and sufficient conditions on the source amplitudes 

In this section we will define the constraints on the active source coefficients B^ n and Bm,n 
by expressing the active source field in terms of near-field and far-field source amplitudes 
and using Graf theorem (|23p . When jxj > |y| the components of can be defined as a sum 
of multipoles at the origin using the first identity in (j23|) 



n=—oo 
oo 

* d = E 4 s v n +(^ X ), 



for Ixl > max(|x m | + a r 



(46) 



where 



M oo 



^ = E E = E E 



M oo 



(47) 



m=l i=— oo 



m=l l=—oo 



If the active field and Vl/^ does not radiate into the far-field, then we must have F^f 1 = 
0, Fn^ = 0, Vn ensuring the necessity of (fl2]) i p. Sufficiency is guaranteed by substituting 
the expressions (fl~2]) i p into an assumed far-field of the form (|4"6j) . 

Next we consider the near-field. Assuming |x m | > a m Mm and using the general form of 
an incident field given by (|6ap . the near-field source amplitudes can be obtained as 



for Ixl < maxflxr 



(48) 



where 



M oo 



^ = E E WW, 



"1 = 1 i = — oo 



M oo 

4 s) = E E 

m=l l=— oo 



(49) 



If the total field is zero in the near-field, then we must have En^ + A^ and + 



ensuring the necessity of (IT2[l ^ /i. Sufficiency is guaranteed by substituting the expressions 
(I12jh a into an assumed near-field of the form (148 p . 
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3.6 Divergence of the active field summation 

The infinite sum expression for the active source fields defined by (|6b|) with source amplitudes 
(jlOap -f fTOcj) is formally valid only in |x — x m | > a m . That is, the expression is not itself valid 
in the domain A in which the sources reside! A valid form could be obtained by using the 
alternative version of Graf's addition theorem in the domain A m associated with the arc 
dC m but the usual form of Graf in the domain A m associated with all other dC n ,n ^ in. 
We would then be assured that the active field is zero everywhere outside C. However if we 
were to do this, the mth source would not be present in the domain A m since the active field 
would be bounded by construction. 

Active cloaking therefore requires that we use the expression (|6bj) with source amplitudes 
(jl0a|) - fjl0c|) for the active field everywhere but we must take a finite number of terms in the 
multipole expansion. That is, we use the source amplitudes that appear in the infinite sum 
as motivation for the choice of source amplitudes that should be chosen in an active field that 
contains only a finite number of multipoles. This ensures a finite (but large) field inside A. 
We should note that this type of difficulty and the fact that it may be used to our advantage 
in the anti-sound context was noted by Kempton [15j . 

With a finite sum for the active field therefore, the integral equation (I20p is not perfectly 
satisfied but instead 



and the field is large (but finite) inside A. Finally we note that a straightforward truncation 
of the active field may not be optimal in terms of cloaking and ensuring a non-radiating field. 
This issue will be considered elsewhere. 

4 Numerical examples 



Figure 2: Plane wave insonification of the cloaking region C generated by four (M = 4) active 
sources placed at the corners of a square. 

The numerical calculations for active source configurations of the type shown in Fig. 
[21 are performed for plane longitudinal and transverse incident waves of a unit amplitude, 
(A p = 1, A s = 0) and (A p = 0, A s = 1), for angles of incidence ip p = ip s = 7°. Variable 
values are taken for the wavenumbers k p and k s , the number of sources M, and the number 
of terms N in summations (|47p and (|49p (the truncation size). The M active sources are 
symmetrically located on a circle of radius b, with 




xeC, 



(50) 




x 



-m 




m = l,M 



(51) 
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where f3 m is the argument of vector x m , and a > b sin 4j . The circular arcs are defined by 



sin 



7T 



7T 



a Sm Mj"M 



m = 1,M. 



(52) 



In all examples, we take a = bs'mjj and consider an elastic medium having a property of 
aluminum with c p = 6427m/s, c s = 3112m/s, p = 2694fc(7/m 3 



4.1 The scattering amplitudes 



|p(p)app| 




Figure 3: Variation of the far-field amplitude coefficients with number of active sources 
(M = 3, 14) for transverse incident waves. In all cases N = 100, ip s = 7°. 

Consider the truncated versions of the infinite sums in eq. (|49p for the farfield amplitudes 
Fn^ and Fn^ , and eq. (|4T|) for the nearfield amplitudes E^ and ■ 



^(p)app ' 



E 



(p)app 



F^ 



(s)app ' 



E, 



{s)app 



M N 



E E < 



m=l l=—N 
M N 

E E a» 

m=l l=-N 



U n-l( X rn), 



Vn G Z. 



Vn G 



(53a) 
(53b) 



The approximate near-field E^ app , En^ app and far-field F^ app , p^ s ' app amplitudes are 
calculated at the incident shear wavenumber k s = 5 varying the number of active sources 
M and the truncation size N. The dependence of the far-field coefficients \F^ app \, \Fn^ app \ 
is illustrated in Fig. [3] for transverse and in Fig. [5] for longitudinal wave incidences. As M 
increases, the far-field coefficients fluctuate at small \n\, and decrease at larger values of |n| 
for both compressional and shear incident waves. 

The variation of the near-field coefficients \An +En p \, \A% -\-En p \ with the number 
of sources M is depicted in Fig. [5] for transverse and in Fig. [6] for longitudinal incident 
waves. For longitudinal wave incidence, the near- field \A ip) + E (p)app \ is less than 10" 4 and 
\An + En^ app \ is less than 10~ 7 . On the contrary, the results are less accurate for transverse 




Figure 4: Variation of the far-field amplitude coefficients with number of active sources 
(M = 3, 14) for longitudinal incident waves. In all cases N = 100, iJj s = 7°. 



waves, as the near-field \A^ + E^ app \ approaches the order of 10 1 and \An + En^ app \ 
reaches a value 10 -4 . 

Figure [7] displays the near-field amplitude coefficients \A$ + En^ app \ and \An + En^ app \ 
as functions of n, the order of Bessel function, for different values of N and M. The accuracy 
of the near-field coefficients improves as N and M increase. Increasing the number of sources 
M allows a decrease in the truncation size N and the order of error. 




Figure 5: Dependence of the near-field amplitude coefficients on n, the order of Bessel func- 
tion, varying the number of active sources [M = 5, 14) for transverse wave incidence. In all 
cases N = 100, = 7°. 
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Figure 6: Variation of the near- field amplitude coefficients with number of active sources 
(M = 5, 14) for longitudinal incident waves. In all cases iV = 100, ip s = 7°. 



4.2 Total displacement field 

4.2.1 Longitudinal plane wave incidence 

First, consider longitudinal plane wave incidence of the form (|3ip . The total displacement 
vector components in Cartesian coordinates are 



\ x i vj \ Q x dx dy ' dy dy dx 
Introducing eq. (j3T|) and (fob)) into (f5"4"]) yields 

M oo 



(54) 



Ux 



I a ir+'h i \\ ■ ■ a ^" f (^p( x x m)) 
B m,n [ cos (fc p (x - x. m ) j - in sm d m — — 



-E E 

m=l n=—c 

+ B$ n I k sin 6» m y n +/ (fc s (x - x m )) + in cos 



V^fc^x - x m )) 



kp | x x r 



+ i cos ipp&i, (55a) 



M oo 



jB$^ n ( sin 6 m V+' (fep(x - x m )) + in cos 6> r 



P m=l n=— oo 



+ ( - K cos 6 m Vn' (k s (x — x m )) + in sin r 



y+(fe p (x-x m )) 



fcp | x x m I 



y+(fc s (x-x m )) 



/up I X x^ 



+ isin-i/'p ^i, 



where 



( ,(x) = arg(x - x,; 



(55b) 
(56) 



4.2.2 Transverse plane wave incidence 

Transverse incident plane waves are of the form (|40p . The total displacement vector compo- 
nents in Cartesian coordinates are 



dx dy dy ' dy dx dx 



(57) 
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Figure 7: The near- field amplitude coefficients as a function of n, the order of Bessel function 
for different values of the truncation size iV in (j53[) generated by different numbers of active 
sources: (a) M = 5 and (b) M = 8, for longitudinal incident waves. 



k s 



Introducing eq. (jlOj) and (l&bj) into (|57j) yields 

M oo 



B m]n \ K 1 cos drnVn' (M X ~ X «0) ~ in sin ^ 



E E 

m=l n=—oo 



+ B^ n ( sin m V+' ( k s (x - x m )) + in cos 8 r 



V+(k p (x-x m )) 



V+(k s (x - x m )) 



k s I x x m I 



kc x x r 



+ i smips^i, (58a) 



M oo 

E E 

m=l n=— oo 



( k 1 sin # m F n + ' (/c p (x - x m )) + in cos 6 r 



V+{k p ( 



/l o X x r 



y n + (fe g (x - x m )) 

T I 

/C s | x x rn i 



i cos vps^i, (58b) 



+ B m,n ( _ cos ^m^'(fcs(x - x m )) + m sin 9 r 
where 6 m is defined by (|56|) . 
4.2.3 Results 

The magnitude of the displacement vector components u x and tt^ are evaluated for ip p = 7° 
for various values of the truncation size N, the number of sources M, and the compressional 
wavenumber k p . Greater accuracy is observed, as expected, with increased TV" and M. How- 
ever, large N and M require longer computation time, and some numerical experimentation 
is necessary to find the smallest values for which the displacement field vanishes to the desired 
degree in the cloaked region. 

The magnitudes of \u x \/k p and \u y \/k p are depicted in Fig. [Hand Fig. [9] for longitudinal 
incidence at different values of N when cloaking devices are active with M = 3, k p = 2. As 
expected, the increase of N is accompanied with the reduction of magnitudes \u x \/k p and 
\u y \/k p in the cloaked region. 

Figure [10] illustrates \u y \/k p for longitudinal incidence with N = 5 changing the values 
of k p and M whilst Fig. [IT] and Fig. [12] show corresponding values of \u x \/k s and \u y \/k s 
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Figure 8: Absolute value of displacement vector component \u x \/k p for N = 5 (a), N = 10 
|(b)[ N = 20 (c) and N = 50 |(d)| when cloaking devices are active with M = 3, k p = 2 for 
longitudinal wave incidence. 



for shear incidence, varying N and M with k p = 2 for the former, and altering the values 
of N and k p with M = 3 for the latter. Comparison of these results shows that at higher 
frequencies, i.e., larger values of k p , greater accuracy is achieved by increasing the number of 
sources M, whereas at lower frequencies the smallest number of sources required, i.e. M = 3, 
produces reasonable cloaking, although this is enhanced with increased values of N. 



5 Conclusions 

The external active acoustic cloaking model of [9] has been generalized to elastodynamics. 
Just as in the former case it is possible to represent the sources in exact terms, although it 
requires that the incident elastic wave field is known in exact form; this is the price paid for 
active control. The control method proposed is based on representing the incident field in 
terms of regular functions (Bessel functions) at each source position, which leads to a linear 
system of equations for the source amplitudes that can be solved in closed form. The linear 
nature of the solution of this essentially inverse problem means that arbitrary incident wave 
motion can be treated by superposition. 

The results presented here provide a first step in the direction of realistic active control 
of elastic waves. Applications to structure borne waves, surface waves, and even geophysical 
waves, are possible. However, as a control problem, many issues remain to be addressed. 
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Figure 9: Absolute value of displacement vector components \u y \/k p for N = 5 (a) , ./V = 10 
|(b)[ N = 20 (c) and N = 30 |(d)| when cloaking devices are active with M = 3, k p = 2 for 
longitudinal wave incidence. 



Not least is the issue of how to balance the goal of silencing one region of space with the 
unavoidable source noise that must be generated in another, larger, region. This quandary 
arises from the fact that the infinite series for the multipole expansion of the m active source 
is divergent inside the domain A m . Exact field cancellation is not achievable in practice; 
it becomes necessary to truncate the series and balance the decrease in cloaking accuracy 
with whatever amplitude level is deemed acceptable in the source region. This is obviously a 
crucial aspect and one that remains to be studied in detail. We have pointed out some simi- 
larities with parallel issues in active noise control, and future studies will examine analogies 
in these topics. One area for consideration is the low frequency end of the spectrum. The 
numerical simulations presented here indicate that a small number of multipoles provide ad- 
equate cancellation at low frequencies. This suggests a natural way to extend ideas based on 
monopoles to more elaborate sources distributions composed of finite numbers of multipoles 
of low order. Hopefully, the present results provide a means to establish realistic strategies 
for practical application. 
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Figure 10: Absolute value of displacement vector component \u y \/k p for k p = 10, M = 3 (a' 
k p = 10, M = 3 |(b)[ k p = 5, M = 7 (c) and A; p = 10, M = 7 (d) when cloaking devices are 



active with N = 5 for longitudinal wave incidence. 



References 

[1] U. Leonhardt. Optical conformal mapping. Science, 312(5781):1777-1780, June 2006. 

[2] J. B. Pendry, D. Schurig, and D. R. Smith. Controlling electromagnetic fields. Science, 
312(5781):1780-1782, 2006. 

[3] S. A. Cummer and D. Schurig. One path to acoustic cloaking. New J. Phys., 9(3):45+, 
2007. 

[4] D. A. Miller. On perfect cloaking. Opt. Express, 14(25): 12457-12466, December 2006. 

[5] F. G. Vasquez, G. W. Milton, and D. Onofrei. Active exterior cloaking for the 2D 
Laplace and Helmholtz equations. Phys. Rev. Lett., 103:073901, 2009. 

[6] F. G. Vasquez, G. W. Milton, and D. Onofrei. Broadband exterior cloaking. Opt. 
Express, 17:14800-14805, 2009. 

[7] F. G. Vasquez, G. W. Milton, D. Onofrei, and P. Seppecher. Transformation elastody- 
namics and active exterior acoustic cloaking. In S. Guenneau and R. Craster, editors, 
Acoustic Metamaterials: Negative Refraction, Imaging, Lensing and Cloaking, pages 1—1. 
Canopus Academic Publishing and Springer SBM, 2012. 



Active elastodynamic cloaking 



November 26, 2012 



21 




n 



I 




I 



-2 -1.5 -1 -0.5 0.5 1 1.5 



-2 -1.5 -1 -0.5 0.5 1 1.5 



(a) N = 5, M = 3 



(b) N = 50,M = 3 



f ... 1^ 

• . it ■ » _ 

• If • » 



-15 I 
2 I 



n 



i 



-2 -1.5 -1 -0.5 0.5 1 1.5 




-1.5 -1 -0.5 0.5 1 1.5 



(c) N = 5, M = 6 



(d) iV = 50, M = 6 



Figure 11: Absolute value of displacement vector component \u x \/k s for N = 5, M = 3 



N = 50, M = 3 [(b)J N = 5, M = 6 [(c)] and TV = 50, M = 6 [(d)] when cloaking devices are 



active with h n = 2,k H = 4.1305 for transverse wave incidence. 



[8] F. G. Vasquez, G. W. Milton, and D. Onofrei. Exterior cloaking with active sources in 
two dimensional acoustics. Wave Motion, 49:515-7524, 2011. 

[9] A. N. Norris, F. A. Amirkulova, and W. J. Parnell. Source amplitudes for active exterior 
cloaking. Inverse Problems, 28 (10) 105002, 2012. 

[10] D. Guicking. On the invention of active noise control by Paul Lueg. J. Acoust. Soc. 
America, 87:2251-2254, 1990. 

[11] P. A. Nelson and S. J. Elliott. Active control of sound. Academic Press, London, 1992. 

[12] R. J. McKinnell. Active vibration isolation by cancelling bending waves. Proc. Roy. Soc. 
A, 421:357-393, 1989. 

[13] C. R. Fuller, S. J. Elliott, and P. A. Nelson. Active control of vibration. Academic Press, 
London, 1996. 

[14] A. David and S.J. Elliott. Numerical studies of actively generated quiet zones. Appl. 
Acoustics, 41:63-79, 1994. 

[15] A. J. Kempton. The ambiguity of acoustic sources - a possibility for active control. J. 
Sound Vibration, 48:475-483, 1976. 



Active elastodynamic cloaking November 26, 2012 22 

[16] G. W. Milton, M. Briane, and J. R. Willis. On cloaking for elasticity and physical 
equations with a transformation invariant form. New J. Phys., 8:248-267, 2006. 

[17] M. Brun, S. Guenneau, and A. B. Movchan. Achieving control of in-plane elastic waves. 
Appl. Phys. Lett, 94:061903, 2009. 

[18] A. N. Norris and A. L. Shuvalov. Elastic cloaking theory. Wave Motion, 49:525-538, 
2011. 

[19] W. J. Parnell. Nonlinear pre-stress for cloaking from antiplane elastic waves. Proc. R. 
Soc. A, 468:563-580, 2012. 

[20] W. J. Parnell, A. N. Norris, and T. Shearer. Employing pre-stress to generate finite 
cloaks for antiplane elastic waves. Appl. Phys. Lett., 100:171907, 2012. 

[21] A. N. Norris and W. J. Parnell. Hyperelastic cloaking theory: Transformation elasticity 
with prestressed solids. Proc. Roy. Soc. A, in press, 2012. 

[22] M. Abramowitz and I. Stegun. Handbook of Mathematical Functions with Formulas, 
Graphs, and Mathematical Tables. Dover, New York, 1974. 

[23] F. Honarvar and A. N. Sinclair. Scattering of an obliquely incident plane wave from a 
circular clad rod. J. Acoust. Soc. Am., 102:41-48, 1997. 



Active elastodynamic cloaking 



November 26, 2012 



23 




-2 -1.5 -1 -0.5 0.5 1 1.5 




-2 -1.5 -1 -0.5 0.5 1 1.5 



(a) N = 5, k p = 2 
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Figure 12: Absolute value of displacement vector components \u y \/k s for N = 5,k p = 2 (a 



N = 5,k p = 5[(b)J A = 20, fc p = 2 [(cJI iV = 20, fc p = 5 p)| A = 50, fc p = 2 [(eJI and 
iV = 50, k p = 5 (f) with M = 3 active sources for transverse wave incidence when cloaking 



devices are active. 
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